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ABSTRACT 

Classic  flutter  flight  testing  involves  the  evaluation  of  a  given  configuration  at  a  stabilized  test  point 
before  clearance  is  given  to  expand  the  envelope  further.  At  each  stabilized  point  flight  test  data  are 
compared  with  computer  simulation  models  to  assess  the  accuracy  of  predicted  flutter  boundaries. 
Because  of  the  time  constraints  associated  with  these  procedures,  the  Air  Force  has  been  seeking  methods 
to  improve  current  flight  test  methods.  This  paper  describes  a  technique  that  provides  a  rapid,  on-line  tool 
for  the  identification  of  aeroservoelastic  systems.  The  technique  involves  the  use  of  discrete  wavelet 
transforms  to  compute  the  impulse  response  (Markov  parameters)  of  the  estimated  system.  This  is  then 
used  in  the  Eigensystem  Realization  Algorithm  (ERA)  method  to  compute  the  discretized  state-space 
matrices.  The  technique  used  herein  includes  metrics  that  are  used  to  assess  the  validity  of  the  identified 
system.  Although  the  method  does  require  that  the  identification  begin  from  stabilized  initial  conditions, 
it  has  been  shown  to  be  relatively  insensitive  to  input  forcing  function.  A  model  of  a  modem  naval  fighter 
aircraft  was  used  to  evaluate  the  capabilities  of  the  identification  method.  The  identification  techniques 
were  evaluated  with  and  without  an  active  oscillation  controller  in  place. 


INTRODUCTION 

Under  certain  common  loading  configurations,  the  F-16C/D  experiences  a  5  Hz  limit  cycle  oscillation 
(LCO)  that  involves  an  interaction  of  asymmetric  wingtip  pitching  and  store  pitching  modes  (Ref.  1). 

This  limited  amplitude  form  of  flutter  is  most  commonly  found  in  the  350  to  450  KCAS  region. 

Although  this  form  of  flutter  does  not  impact  structural  integrity,  it  can  produce  significant  ride  quality 
problems  for  the  pilot.  Lockheed  Martin  Tactical  Aircraft  Systems  was  contracted  by  the  Air  Force  to 
develop  and  demonstrate  a  non-adaptive  active  flutter  suppression  system  to  reduce  or  eliminate  this  LCO 
phenomenon.  The  resulting  control  laws  use  the  flaperons  to  counteract  the  motion  of  the  5  Hz  LCO. 
Although  the  system  was  not  successful  at  all  flight  conditions  and  loading  configurations  tested  in  flight, 
it  did  demonstrate  a  means  for  allowing  an  expanded  flight  envelope  on  an  operational  aircraft. 

According  to  Ref.  1  specific  problems  encountered  in  the  flight  test  program  included  the  requirement  for 
a  very  detailed  and  complex  test  procedure  to  optimize  AFS  control  gains.  In  addition,  poor  initial  gain 
and  phase  predictions  for  the  AFS  were  not  well  suited  to  suppress  the  LCO  and  the  required  gains 
changed  quickly  as  the  aircraft  accelerated  above  Mach  0.9  to  0.95. 

This  program  and  other  related  work  sponsored  by  the  Air  Force  has  emphasized  a  need  for  improved 
aeroservoelastic  prediction  models  and  flutter  flight  test  techniques.  In  response  to  the  need  for  improved 
models,  significant  research  and  development  has  been  conducted  by  the  University  of  Colorado  at 
Boulder  (Refs.  2  and  3).  The  approach  taken  is  to  use  nonlinear  methods  that  better  capture  the  unsteady 
flow  characteristics  in  the  transonic  region  that  features  mixed  subsonic-supersonic  flow  patterns.  The 
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technique  has  been  applied  to  stabilized,  accelerated,  and  increased  angle-of  attack  configurations. 
Systems  Technology,  Inc.  has  been  contracted  by  the  Air  Force  Flight  Test  Center  to  develop  improved 
flutter  flight  test  techniques  for  aircraft  with  and  without  a  flutter  suppression  system.  This  work  has  led 
to  the  Wavelet-Eigensystem  Realization  Algorithm  system  identification  method,  among  others,  that  is 
the  focus  of  this  paper. 


THE  WAVELET-ERA  METHOD 

A  methodology  has  been  developed  that  estimates  the  frequency  response  of  a  given  system  using  an 
arbitrary  input  and  the  system’s  response.  This  technique  uses  the  discrete  wavelet  transform  as  given  in 
previous  progress  reports.  However,  if  it  can  be  assumed  that  the  system  is  initially  at  rest  (i.e.,  the 
aircraft  is  trimmed),  then  the  technique  is  greatly  simplified.  When  the  initial  conditions  are  zero,  the 
system  response  consists  entirely  of  a  forced  response  from  the  time  the  input  is  first  applied  until  the  end 
of  the  data  measurement.  Thus  it  is  not  necessary  to  estimate  what  the  system  dynamics  are,  which  would 
have  been  required  for  an  estimate  of  the  free  response. 

The  discrete  wavelet  transform  of  the  impulse  response  is  computed  using  the  forcing  input  and  system 
response,  and  the  inverse  wavelet  transform  is  then  computed  to  produce  the  time-domain  impulse 
response  (also  known  as  the  Markov  parameters).  Finally,  the  Eigensystem  Realization  Algorithm 
(ERA)  method  is  used  to  compute  the  discretized  state-space  matrices,  from  which  is  generated  the 
magnitude  and  phase  of  the  system  response. 

The  experimental  identification  of  physical  vibration  modes  and  mode  shapes  can  be  done  with  impulse 
response  functions  that  are  extracted  from  measured  vibration  records.  In  order  to  extract  these  impulse 
response  functions,  the  Fast  Fourier  Transform  (FFT)  has  often  been  used  in  conjunction  with  repeated 
data  filtering  and  windowing  (Ref.  4).  These  techniques  require  input  signals  that  are  rich  enough  to 
excite  structural  frequencies  of  interest  -  this  does  not  appear  to  be  the  case  with  wavelet  analysis. 
Moreover,  with  its  more  intuitive  decomposition  of  data,  wavelets  analysis  allows  identification  of  time- 
varying  system  parameters. 

The  system  realization,  based  on  a  first-order  state  space  model,  can  be  represented  as: 

x( t)  =  Axft)  +  Bu(t) 
y  =  Cx(t) + Du(t) 


In  the  time  domain,  the  solution  at  time  tn  to  Equation  1  can  be  expressed  as: 

y(U-  J '”Ji(tn-T)u(r)dT 


(1) 


(2) 

where  h(t)  is  the  temporal  impulse  response  function.  This  convolution  formula  can  be  expressed  in 
matrix  form  by: 

Y  -  hU 


(3) 

In  Equation  3,  Y  represents  the  output  matrix,  U  is  the  input  matrix,  and  h  is  the  time-domain  impulse 
response  matrix.  Once  h  is  known  the  output  of  our  system  can  be  determined  for  any  arbitrary  input.  An 
accurate  extraction  of  h(t),  often  referred  as  the  Markov  parameters,  will  identify  a  system  (Ref.  5)  since 
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h(t)  =  CeA'B 

(4) 

One  could  try  to  obtain  h(t)  directly  by  applying  a  unit  impulse  input  u(t).  In  this  case,  y(t)  would 
theoretically  be  identical  to  our  desired  h(t).  However,  this  is  practically  unfeasible.  Looking  at  Equation 
3,  if  the  input  signals  are  not  rich  in  frequency  content,  or  the  sampling  size  is  too  large,  U  can  become  ill- 
conditioned  and  the  corresponding  impulse  response  function  cannot  be  computed  accurately.  Moreover, 
the  computations  associated  with  h(t)  can  be  intensive.  For  these  reasons,  FFT-based  extraction  of  the 
impulse  response  function  is  widely  used  (Ref.  4)  since  it  has  high  computational  efficiency  and 
reasonable  robustness  provided  the  input  data  is  rich  in  frequency  content.  However,  if  the  input  load  u 
consists  of  only  a  single  or  a  few  frequencies,  the  temporal  impulse  response  data  tends  to  be  erratic  and 
badly  conditioned.  Wavelet  analysis  appears  to  offer  a  more  robust  alternative. 

COMPUTATIONAL  TECHNIQUE 


Wavelet  Theory 

The  wavelet  transform  allows  any  arbitrary  signal  f(x)  to  be  decomposed  into  an  infinite  summation  of 
wavelets  at  different  scales  according  to  the  expansion  (Ref.  6): 

m=  £  £cMw(vx-k> 

j=- co  *=- oo 

(5) 

The  W(2lx-k)  functions  are  the  wavelets  and  provide  a  local  basis  along  the  time  axis.  To  see  how 
wavelets  can  be  generated  from  the  so-called  dilation  equation  and  their  relation  to  the  scaling  functions 
(p(x),  the  reader  is  directed  to  Ref.  7.  Because  of  the  way  in  which  the  wavelets  are  defined,  when  j  is 
negative,  the  wavelets  W(2’x-k)  can  be  expressed  as  a  sum  of  terms  like  (p(x-k).  The  corresponding 
wavelet  transform  can  then  be  expressed  as  follows: 

f(x)=  YJc<pk(p(x-k)+Yj  ±chkW(Vx-k) 

jfc=-oo  y=o  *=-oo 

(6) 

In  this  approach,  the  Daubechies  coefficients  (Refs.  8  and  9)  are  used  to  generate  these  wavelets.  Note 
that  db2,  which  is  used  in  this  analysis,  spans  a  bit  less  than  3  units  over  the  x-axis.  Daubechies’  family  of 
wavelets  satisfies  two  crucial  requirements:  orthogonality  of  local  basis  functions,  and  second  or  higher- 
order  accuracy,  depending  on  the  dilation  expression  adopted  (Ref.  10). 

In  order  to  use  the  wavelet  transformation  in  practical  applications,  a  way  to  define  a  discrete  version  of 
the  transformation  shown  in  Equation  7  must  be  found.  For  this  purpose,  the  range  of  the  independent 
variable  x  is  limited  to  the  unit  interval  [0,1),  and  a  wavelet  series  expansion  performed  over  that  interval 
(Ref.  1 1).  A  complication  arises  since  some  of  the  wavelets  W(2x-k)  overlap  the  edges  of  the  interval. 

For  this  reason,  it  is  convenient  to  assume  that  f(x)  is  one  period  of  a  periodic  signal  exactly  repeated  in 
the  adjacent  unit  intervals.  Note  that  with  the  D4  wavelet  family  W(x)  spans  almost  three  units.  Over  the 
interval  [0,1)  there  are  contributions  from  three  bases:  from  the  first  third  of  W(x),  from  the  middle  third 
of  W(x+1),  and  from  the  last  third  of  W(x+2).  This  is  equivalent  to  W(x)  being  wrapped  around  the  unit 
interval.  Thus  the  wavelet  expansion  of  f(x)  in  the  [0,1)  interval  can  be  written  as: 
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(W(2x) 
W(2x-l)\ 


+  [«4 


a. 


W(  Ax) 
W(4x-\) 
W(4x~2) 
W(4x-3) 


+  alJ+kW(  Vx-k)  +  ... 

(7) 

The  coefficients  ai,  a 2,  a3,  84,  ...  give  the  amplitudes  of  the  contributing  wavelets  (after  wrapping)  to  one 
cycle  of  the  periodic  function  in  the  [0,1)  interval.  Because  of  orthogonality  conditions,  the  general 
wavelet  transform  coefficients  can  be  found  by 


a11+k=2J\f(x)W(Vx-k)dx, 

0 


1 

«o  =  \f(x)(p(x)dx 
0 


(8) 

The  discrete  wavelet  transform  (DWT)  is  an  algorithm  for  computing  Equations  7  and  8  when  f(x)  is 
sampled  at  equally  spaced  intervals  over  [0,1).  As  stated  before,  it  is  assumed  here  that  f(x)  is  one  period 
of  a  periodic  signal  and  that  the  scaling  and  wavelet  functions  wrap  around  the  interval  [0,1).  A 
remarkable  feature  of  the  DWT  algorithm  is  that  there  is  no  need  to  compute  cp(x)  or  W(2x-k)  explicitly. 

A  MATLAB  implementation  of  this  DWT  algorithm,  found  in  Appendix  7  of  Ref.  1 1 ,  was  discovered  by 
Mallat,  and  is  often  referred  by  as  the  Mallat ’s  pyramid  algorithm. 

DWT-Based  Extraction  of  the  Markov  Parameters 

The  DWT  method  starts  with  the  convolution  integral, 

t„  1 

y(t„)=  \h(tn  -T)u(r)dz  =  \h(0)u(tn  -6)d6 

-00  0 


(9) 

Note  that  h(6)  is  the  temporal  impulse  response  function.  The  impulse  response  function  is  expanded  in 
terms  of  the  wavelet  basis  functions  for  the  entire  time  interval,  0  <  9  <  1 ; 

m =h2WT + Z  Z  hDZw(2i  9+k) 

j  * 


(10) 

where  the  hhWi  terms  are  the  expansion  coefficients.  For  the  DWT  characterization  of  u(t„-9),  first  u(9)  is 
reversed  in  time  to  obtain  u(-9),  then  it  is  shifted  toward  the  positive  time  axis  by  an  amount  t„.  Following 
this,  u(t„-9)  can  be  expressed  as 

u(tn  -0)  =  u7T  +  IX  uJIFW  e  +  ^ 


01) 

Substituting  Equations  10  and  1 1  into  Equation  9,  and  making  use  of  the  orthogonality  conditions,  yields 
the  following  formula  (Ref.  1 1): 
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y{tn)  =  hDWTuDm(tn) 


For  the  entire  data  sample,  the  inputs  and  outputs  are  arranged  in  the  form 

y  —  fjDWTjjDWT 

Solving  for  hDm  (Ref.  9), 

ft dwt  _  y  jjDWTt  (jjdwt  ijDwrTyi 


(12) 


(13) 


(14) 

The  desired  temporal  impulse  response  data  is  finally  obtained  by  taking  the  inverse  wavelet  transform  of 
hDWT.  (See  Ref.  1 1  for  details  on  how  to  implement  the  inverse  wavelet  transform.) 

h(t)  =  IDWT(hDWT ) 


(15) 


ANALYSIS  EXAMPLE 


Aircraft  Model  Description 

During  developmental  flight  testing  of  a  modem  fighter  aircraft,  an  unacceptable  5.0  to  6.0  Hz  LCO  was 
experienced  at  high  speed,  low  altitude  flight  with  certain  store  configurations  (Ref.  12).  The  approaches 
investigated  to  suppress  the  limit  cycle  included  variations  in  store  orientation,  wing  control  surface 
biasing,  and  active  oscillation  suppression.  Changes  in  store  orientation  and  biasing  techniques  did  not 
sufficiently  reduce  the  problem,  so  various  active  suppression  techniques  were  evaluated.  The  successful 
active  oscillation  control  approach  used  a  passive  band  pass  filter  to  command  aileron  deflections  with 
gain  and  phase  set  by  flight  control  system  lateral  acceleration  input  signals.  This  system  was 
incorporated  into  the  production  aircraft  using  existing  flight  control  system  components. 

Under  a  subcontract  to  STI,  The  Boeing  Company  provided  aeroservoelastic  models  of  this  fighter 
aircraft  (see  also  Ref.  13).  The  aircraft  was  configured  with  full  330-gallon  tanks  on  the  inboard  pylons 
and  MK-84s  (i.e.,  2,000  lb  bombs)  on  the  outboard  pylons,  a  configuration  that  produced  LCO  in  flight. 
The  models  were  formulated  at  twelve  flight  conditions,  below,  at,  and  above  the  aircraft’s  flutter 
boundaiy,  as  identified  in  Figure  1 .  The  lateral-directional  model  was  provided  in  the  form  of  ABCD 
state  space  matrices  for  the  rigid  body  airframe,  aeroservoelastic  (ASE)  dynamics,  rigid  body  control 
system,  ASE  control,  plant  input  time  delay,  zero  order  hold,  and  sensors.  Inputs  to  the  model  are  the 
stick  and  rudder  positions,  and  outputs  consist  of  the  rigid  body  dynamics,  flexible  dynamics,  and 
accelerator  data  at  various  locations  on  the  aircraft.  A  block  diagram  of  the  model  is  shown  in  Figure  2. 
The  ASE  dynamics  and  the  active  oscillation  controller  (AOC)  can  be  turned  on  or  off.  A  block  diagram 
and  Bode  frequency  response  of  the  band  pass  filter  implementation  of  the  AOC  is  shown  in  Figure  3.  As 
described  in  Ref.  13  the  controller  is  activated  as  a  function  of  Mach  number  and  altitude,  and  only  for  the 
specific  set  of  stores  where  LCO  is  known  to  occur. 
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Figure  1.  Aeroservoelastic  Model  Flight  Conditions 
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Figure  2.  Modern  Fighter  Lateral  Axis  Rigid  Body  and  ASE  Dynamics,  with  Rigid  and  AOC  Control 
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a)  Filter  Block  Diagram:  Input  -  ny  Aoc  (g)»  Output  -  8a  (deg) 


Frequency  (rad/sec) 
b)  AOC  Frequency  Response 

Figure  3.  Active  Oscillation  Controller 
Assessment  of  WERA  Identification 

Figure  4a  compares  the  frequency  responses  of  the  actual  and  estimated  aeroservoelastic  dynamics  for  the 
Mach  0.9,  10,000  ft  altitude  flight  condition,  showing  excellent  agreement  in  both  magnitude  and  phase. 
The  discrete  wavelet  transform  is  used  to  compute  the  impulse  response  (Markov  parameters)  of  the 
estimated  system,  which  is  used  in  the  Eigensystem  Realization  Algorithm  (ERA)  method  (Ref.  14)  to 
compute  the  discretized  state-space  matrices.  Accurate  system  identification  thus  rests  in  large  part  on  the 
quality  of  the  impulse  response  estimate.  In  Figure  4b  the  actual  and  estimated  impulse  responses  are 
practically  indistinguishable  from  one  another;  however,  there  may  be  cases  where  the  estimated  impulse 
response  degrades  significantly  toward  the  tail  of  the  impulse  data-window  (due  to  noise,  precision  error, 
etc.).  This  degradation  typically  is  seen  as  a  rapid  oscillatory  divergence  in  amplitude.  Normally  the 
impulse  response  sees  a  maximum  peaking  shortly  after  start-up,  followed  by  decaying  oscillations  that 
taper  towards  zero  or  resonate  with  low-amplitude.  In  order  to  ensure  proper  choice  of  impulse  length, 
the  impulse  response  is  divided  into  time  bins  following  the  maximum  start-up  peak,  and  the  standard 
deviation  of  these  bins  is  computed.  If  the  standard  deviation  of  a  bin  is  greater  than  the  one  preceding  it, 
the  impulse  response  is  terminated  at  the  preceding  bin.  This  procedure  systematizes  window-length 
selection,  and  ensures  maximum  use  of  available  data. 
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b)  Impulse  Response 

Figure  4.  Model  Estimates  and  Comparisons  for  Mach  =  0.9,  Altitude  =  10,000  ft  Flight  Condition 

In  Equation  16  the  ABCD  matrices  identified  with  the  ERA  method  are  transformed  to  modal  coordinates, 
and  Equation  17  computes  the  modal  time  history  for  the  ;th  mode.  The  modal  pulse  response  used  in  the 
identification  is  computed  with  Equation  18,  where  RES1  are  the  singular  value  decomposition  matrices  of 
the  Hankel  matrix  whose  elements  are  the  impulse  responses.  The  dot  product  between  Equations  17  and 
1 8  is  defined  as  the  Modal  Amplitude  Coherence  (MAC),  shown  in  Equation  19.  The  MAC  is  an 
important  metric  of  confidence  that  compares  the  modal  time  history  computed  two  different  ways  -  if  the 
two  methods  yield  a  close  match  for  a  particular  mode,  there  is  high  likelihood  that  the  mode  has  been 
correctly  identified.  Figure  5  shows  the  MACs  for  this  example  to  be  essentially  unity  at  all  modes, 
which  is  consistent  with  the  excellent  fits  seen  in  Figure  4a  and  b. 

A  =  'PA'F  '  A=  A 

m 

Bm=^'B  Cm=C¥  Dm  =  D 
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q, 


=  [4  X,b,  A,2b,  -  A‘-% 


Q  =  r~%  Sf 


MACi= 7 


m 

q,q  * 

(16) 

(17) 

(18) 


(19) 


Modal  Amplitude  Coherence 


50  100 

Modal  Damping 


Figure  5.  Mode  Amplitude  Coherence  and  Damping  for  the  Figure  4  Model  Estimates 

A  metric  that  indicates  the  relevance  of  the  mode  identified  is  given  in  Equation  20,  called  the  Mode 
Singular  Value  (MSV).  This  is  a  measure  of  the  contribution  of  each  identified  mode  to  the  identified 
pulse  response  history  -  if  a  pole  is  nearly  cancelled  by  a  zero  one  would  expect  the  MSV  at  that  mode  to 
be  very  small.  While  it  is  not  feasible  for  an  identification  algorithm  to  compute  all  modes  of  a  system,  a 
system  has  effectively  been  estimated  when  actual  modes  corresponding  to  high  MSVs  are  closely  aligned 
with  estimated  modes  that  exhibit  similarly  high  MSVs.  The  algorithm  searches  for  the  closest  estimated 
mode  to  every  actual  mode  (using  natural  frequency  and  damping  as  distances).  The  absolute  value  of  the 
MSV  difference  between  the  actual  and  estimate  is  then  computed,  shown  in  Figure  3d  where  for  clarity 
only  damping  values  less  than  0. 1  are  shown.  If  more  than  one  actual  mode  is  in  close  proximity  to  an 
estimated  mode,  the  actual  mode  with  the  highest  MSV  is  differenced  with  the  estimate.  The  MSVs  of 
the  remaining  actual  modes  that  have  not  been  paired  with  estimated  modes  are  simply  plotted  by 
themselves.  A  well-estimated  system  would  thus  exhibit  MSV  differences  and  stand-alone  values  that  are 
all  approximately  the  same  order  of  magnitude,  indicating  that  the  actual  modes  not  identified  were 
comparatively  negligible.  This  in  fact  is  the  case  in  Figure  3d. 
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MSVt  =  ^|c,|  (1 + |i,-|+ |ij2 1 + •  •  • + 1^/_2|)|*, 


(20) 


Figure  6.  Mode  Singular  Value  Differences  for  the  Figure  4  Model  Estimates 

The  MSV  metric  defined  above  is  for  the  identified,  digital  system.  We  have  explored  the  use  of 
equivalent  metrics  for  analog  systems.  This  is  important  because  our  model  is  analog,  and  high  order 
(greater  than  100),  and  numerically  it  is  better  to  compute  the  MSV  using  the  analog  model.  Define  the 
following  analog  system  and  the  ZOH-equivalent  digital  system: 

f  x  =  Ax  +  Bu 
Analog  system:  { 

[y  =  Cx  +  Du 


Eigenstructure:  A  =  XAX~'  ,  X  =  [x]---x„]  ,  A=diag(/?,  •••An)  ,  X~T  =  [y, 


Residue  of  mode  A,  =  a,  +  jco:  :  Rt  =  Cxty]  B 


analog  MSV,  =  prcr^  (R,) 
o', 


{X(fC  |)  — —  ^  |  ^  ». 

,  where  <t>=eAT  ,  F  =  f  eAr Bdr 
y  =  Cx+Du  * 


Eigenstructure:  <I>  =  XeKT X  1  ,  X,A,X  1  =  same  as  for  analog  system 


Residue  of  mode  r,  =ev  :  S,=  Cx,yJ r 


(21) 

(22) 

(23) 

(24) 

(25) 

(26) 
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The  new  result  is  that  the  residues  and  modal  singular  values  are  related  by: 


s. 


digital  MSV,  =  analog  MSV,  x 


J~z\ 

4 

> 

1 

JM 

(27) 
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(29) 


Beyond  the  Flutter  Boundary  Analysis 

Figure  7  shows  results  using  Flight  Condition  6  with  the  Oscillation  Controller  activated.  When  the 
Oscillation  Controller  is  deactivated,  an  unstable  pole  is  identified  at  approximately  34  rad/sec,  which 
agrees  with  the  model  for  this  flight  condition  (Figure  8).  Note  the  divergence  of  amplitude  in  the  Figure 
8  time  histories. 

CONCLUSIONS 

A  method  for  rapid,  on-line  identification  of  aeroservoelastic  systems  has  been  developed  that  employs 
discrete  wavelet  transforms.  The  wavelet  transforms  compute  the  system  impulse  response,  which  is  then 
used  in  the  Eigensystem  Realization  Algorithm  (ERA)  method  to  generate  the  discretized  state-space 
matrices.  The  method  was  demonstrated  using  an  aeroservoelastic  model  of  a  modern  military  fighter. 
Modal  Amplitude  Coherence  and  Mode  Singular  Value  metrics  were  defined.  Excellent  agreement  was 
observed  between  the  aeroservoelastic  model  and  the  estimated  system. 
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Figure  7.  Model  Estimates  and  Comparisons  for  Mach  =  0.9,  Altitude  =  5,000  ft  Flight  Condition  AOC  On 
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c)  Mode  Amplitude  Coherence  and  Damping 

Figure  8.  Model  Estimates  and  Comparisons  for  Mach  =  0.9,  Altitude  =  5,000  ft  Flight  Condition  AOC  Off 
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